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Abstract. In this paper, we prove almost global existence of solutions to certain 
quasilinear wave equations with quadratic nonlinearities in infinite homogeneous 
waveguides with Neumann boundary conditions. We use a Galerkin method to ex- 
pand the Laplacian of the compact base in terms of its eigenfunctions. For those terms 
corresponding to zero modes, we obtain decay using analogs of estimates of Klain- 
erman and Sideris. For the nonzero modes, estimates for Klein-Gordon equations, 
which provide better decay, are available. 



1. Introduction. 

The purpose of this paper is to give a simple proof of almost global existence for 
quasilinear Neumann wave equations on infinite homogeneous waveguides. This work 
expands on results of Metcalfe, Sogge, and Stewart ^D] an d of Lesky and Racke [§]. Here, 
as in these works, the key step will be the use an eigenfunction expansion in the compact 
base and to use estimates for Klein-Gordon equations for those terms corresponding to 
nonzero modes. For the zero modes, we will use decay estimates that are analogous to 
those of Klainerman and Sideris [Hj- In order to prove long time existence, we couple 
these decay estimates with energy estimates. 

Let us describe our initial-boundary value problem more precisely. We will be studying 
nonlinear wave equations in infinite homogeneous waveguides, R 3 x ft where C M. d 
denotes a nonempty, bounded domain with smooth boundary. We examine equations of 
the form 



(1.1) 



Uu = Q(d t>x u, dl x u), (t, x, y) e R+ x M 3 x O, 
d v u(t,x, - )\ d n = 0, 

u{0, x, y) = f(x, y), d t u(0, x, y) = g(x, y), 



where □ = &l — (A + An) is the d'Alembertian on the waveguide. Here, <9„ denotes the 
normal derivative on d£l. Moreover, 



A = A R3 = d2 / dx . 
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is the Laplacian on R 3 and 

3=1 

is the Neumann Laplacian. 

The nonlinearity Q is quadratic in its arguments and is afhne linear in d^ x u (i.e. is 
quasilinear) . We may expand Q as follows 

(1.2) Q(d t , x u,df tX u) = ^2 Af k diudjdku + R(dt, x u,dt,xu) 

0<j,k,l<3 

where the Aj k are real constants and R is a constant coefficient, quadratic form. Here 
and throughout, we set xq = t and do — dt when convenient. 

In order to solve l|l.lf> , the data must be assumed to satisfy the relevant compatibility 
condition. Let JkU = {d xy u : < \a\ < k} denote the collection of all spatial derivatives 
of u of order up to k (using local coordinates in a small tubular neighborhood of dSl). If u 
is a formal H N solution for some fixed N, then we can write d^d v u(0, ■ ) = ^k(Jk+if> Jkg) 
where the ^k are called compatibility functions and depend on Q, Jk+if and Jkg- The 
compatibility conditions for (/, g) E H x H N ^ simply require that the ^>k vanish on 
R+ x R 3 x dfl when < k < N - 2. Moreover, we say that (f,g) E C°° satisfy the 
compatibility conditions to infinite order if this condition holds for all N. 

We assume that the initial data have compact support and are small in norm. That 
is, we assume that there is a fixed constant B > so that 

(1.3) f(x,y)=g(x,y)=0, \x\>B. 
Moreover, we assume that 

(1-4) ||/||j?JV(R3xn) + \\9\\h n - 1 («. ; > x q) < £, 

where 

H/llff«(R 3 xSl) = X! \\dx,yf\\L2(R3xn)- 
\a\<N 

Under these assumptions, we can prove the following almost global existence result. 

Theorem 1.1. Assume that the Cauchy data (/, g) E C°°(R 3 x fi) satisfy and (|1.4|) 
as well as the compatibility conditions to infinite order. Then there are constants N , k 
and £q > so that if e < Eq and N is sufficiently large in (|1.4fl , then (|l.lf> has a unique 
solution u E C°°([0, T £ ) xM 3 x!i) where 

(1.5) T e = cxp(K/e). 

Notice that the lifespan (|1.5|l is sharp. Indeed, if one takes f(x,y) = f(x), g(x,y) = 
g{x) independent of y, then solutions of (|1.1|) are equivalent to solutions to Ou = 
Q(dt. x u, Of x u) in Minkowski space. By the classical counterexamples of John (see, e.g., 
@J0>]) of the form Uu = (dtu) 2 , it is seen that (|1.5fl cannot be improved. See, also, Sideris 

EH- 

It is for technical reasons that we are only able to handle nonlinearities that do not 
depend on the derivatives d y . However, when the compact base is one dimensional, we 
are able to obtain an optimal result. In order to use energy methods, we must now 
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(1.6) 



assume a nonlinear compatibility condition as in |10|. Here, we will be studying the 
initial-boundary value problem 

'Uu = Q(du,d 2 u), (t,x,y) el + x R 3 x [a, b], 
d y u(t,x,a) = 0, d y u(t,x,b) = Q 
u(0, x, y) = f(x, y), d t u(0, x, y) = g{x, y), 

where a < b are fixed constants and d — dt, x ,y denotes the full space-time gradient. 
Expanding Q as above, we have that 

(1.7) Q{du,d 2 u)= ^2 B^diudjdku + Ridu^u), 

0<j,k,l<4 

where R is again a constant coefficient quadratic form. Here, for convenience, we set 
X4 = y and — d y . In order to use energy methods with Neumann boundary conditions, 
we must assume the following nonlinear Neumann compatibility condition 

(1.8) ]T BfO'ij", o. if(6,t,n)ex, 

0<j,k,l<4 

where 

(1.9) X = {(0,Z,t)) ■ e=(o,o,o, v (y)),^e = o,r l -e = o,yGdn}. 

I.e., when 6 is normal to i? 1+3 x dfl and £ and rj are orthogonal to 9, l|1.8l) must hold. 
As this condition automatically holds when the quasilinear terms only involve djdkU, 
< j, k < 3, such an assumption was unnecessary when studying (|1 . 1|) . This is the 
natural condition that is required in order to obtain energy estimates for quasilinear 
equation with Neumann boundary conditions. 

Under these assumptions, we can then prove 

Theorem 1.2. Assume that the data (f,g) G K 3 x [a,b] satisfy (|1.3|l and ljl.4|l as well 
as the compatibility conditions to infinite order. Moreover, assume that Q satisfies (|1.8(l . 
Then, there are constants N , n, and Eq so that if e < Sq and N is sufficiently large in 
(|1.4I) . then (|1.6[) has a unique solution u G C°°([0,T E ) x M 3 x [a, b]). 

Results similar to Theorem 11.11 and Theorem 11.21 were first obtained by Lesky and 
Racke in K+ x R™ x f2 when n > 5. These results were improved in Metcalfe, Sogge, 
and Stewart There, for Dirichlet boundary conditions, global existence was shown 
in M+ xl"x!l where Q C R d for any n > 3. By using an eigenvalue expansion, one can 
reference decay estimates for the Klein-Gordon equation which provide 0(tr n l 2 ) decay 
as opposed to the standard O^ - ^"" 1 )/ 2 ) decay that is available for the wave equation. 

When Neumann boundary conditions are assumed, the proofs are more delicate. First 
of all, in order to use energy methods a condition of the form Ijl.SJI must be assumed. 
This is then sufficient to prove global existence for Klein-Gordon equations. For the wave 
equation, as there is a zero mode, the estimates that are available do not provide the 
same decay. In |TJ]], the techniques of [7] were applied to obtain existence for semilinear 
wave equations whose quadratic nonlinearities only depend on d x u. 

Theorem 11.11 and Theorem 11.21 are improvements to this last result of |10| . Here, we 
instead use the eigenfunction expansion. For those terms corresponding to zero modes, 
we use estimates analogous to those of |5] for the wave equation. For the nonzero modes, 
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we have the better estimates for Klein-Gordon equations. As the estimates of [5] only 
provide decay for the gradient c^u, we are lead to the restriction on the nonlinearity 
in Theorem ll.il Theorem II . 21 follows from the simple observation that d y u = d v u when 
the compact base is one dimensional. Thus, in this case, d y u solves a wave equation 
with Dirichlet boundary conditions. Hence, better decay estimates are available for these 
terms in this special case. 

This paper is organized as follows. In the next section, we gather some weighted 
Sobolev estimates and estimates for the wave equation that are analogs of those in [5| . In 
the following section, we show that our necessary decay estimate for the wave equation 
follows from those estimates in the previous section. We also provide the decay estimate 
that we shall use for solutions to Klein-Gordon equations. This is essentially from 
In the fourth section, we extend the previous estimates (which were on M + x K 3 ) to 
estimates on the waveguide. To do so, we use a Galerkin expansion in the y variable 
and use the inequalities from the previous section. In the fifth section, we examine the 
energy estimates which we shall require. These are rather standard, except that we need 
to make a natural nonlinear compatibility assumption on the quasilinear terms. Such 
an assumption is necessitated by the Neumann boundary conditions. Finally, in the last 
section, we prove Theorem 11.11 and Theorem ll.2l 

It is a pleasure to thank C. Sogge for his collaboration ^U] that preceded this paper 
and for helpful suggestions concerning the current study. 

2. Sobolev estimates and decay estimates for wave equations in R + x K 3 . 

In this section, we present the estimates that provide the necessary decay for wave 
equations in R + x R 3 . These will be used to estimate the terms in the expansion of u 
with zero mode. Here, we let 

(2.1) {T} = {dt,d x ,Q jk : 0<j<k<3} 
where 

= Xidj — Xjdi, 1 < i < j < 3 

and 

^ofc = xkd t +tdk-, 1 < k < 3. 

We will also denote 

z = {d t ,d x ,n jk : i<i<*<3}. 

The first lemma that we shall require is the following weighted Sobolev inequality. 
Lemma 2.1. Let w S (M+ x K 3 ). Then, 

(2.2) <r) 1/2 H^)|< E 

\a\<l 

and 

(2.3) {r)\dT a w{t,x)\< E P*0u;(V)ll2. 

\l3\<\a\+2 
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Here and throughout, we use the notation A < B to indicate that A < CB for some 
positive, unspecified constant C. We are also using (x) — (r) = yvL-j- \x\ 2 . The first 
of these estimates essentially appears in Sideris See also Hidano [2]. The second 

estimate is from Klainerman and Sideris [Hj- As a corollary of (|2.2J) . we obtain 

Corollary 2.2. Let w S C °°(R + x R 3 ). Then, 

(2.4) (r^-^isrM^)! < E H r/V (*>-)ll2+ E ll(i-r)9 2 r^(t, - > || 2 - 

l/3|<M + l l/3|<M + l 

This estimate also appeared in [2] and follows by simply applying 1)2. 2(1 to (t—r)dT a w(t, x) 

In order to estimate the last term in (|2.4|l . we will require estimates that are analogs 
of those due to Klainerman and Sideris 0. . In |Hj, the necessary estimates were restricted 
to only make use of the vector fields Z and the scaling vector field L = tdt + rd r . In 
the current setting, it is more convenient to use estimates that only involve T. Thus, we 
prove the following variant of the estimates of [8]. 

Lemma 2.3. Let w e C 2 (R+ x R 3 ). Then, 



(2.5) (t-r)\Aw(t,x)\ < E \dT a w(t,x)\ + (t + r)\n w (t,x)l 

\a\<l 

(2.6) (t-r)\d?w(t,x)\ < E \dT a w(t,x)\ +t\D w {t,x)\, 

\a\<l 

(2.7) (t-r)\V x d t w(t,x)\< E \dT a w(t,x)\+t\a w {t,x)\. 

M<i 

Here and throughout the current section and the next, we are abusing notation and 
using U = df — Ar3 to denote the d'Alembertian on Minkowski space. 

The proof is based in part on the following lemma, which also appeared in [S]. 
Lemma 2.4. Let w E C 2 (R 3 ). Then, 



(2.8) \Aw(x) -&?w(x)\ < - V \V x Z a w(x)\. 



r 



\a\<l 



This estimate follows immediately from the expansion of the Laplacian into its radial 
and angular parts 

a = d 2 r + ^-d r + 

where il — — A V x and A is the usual vector cross product on R 3 , as well as the fact 
that 

\d r w( •)! < |Vxtu( 01, -Vw( ■ ) < |V«>( 01- 
r 

Proof of Lemma WSi The proof will be broken into the cases that \x\ < t/2 and \x\ > t/2. 

Case 1 (\x\ < t/2): Our main technique here will be the use of substitutions involving 
f^oi = tdi + Xidt, as well as the fact that the commutator of any T with d only involves 
d. 
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In order to obtain (|2.5p . we begin with 



fl 0l diW = ^ (tdi + Xid t )diW 

<i<3 l<i<3 

= tAw + Xididtw 

l<i<3 

X T 2 

= tAw + ^ -j-^oidtw - —dl 

l<j<3 

2 



9 9 

t— — )Aw+ \_. —^oidt.w -Ow. 

l<i<3 

2 

If we solve for the term involving A and use the fact that t — > t — r in this region, 
we see that 

(t-r)\Aw\< \tt idiw\+ Y \tooidtw\+tpw\, 

l<i<3 l<i<3 

from which (|2.5(l follows. 

We use similar arguments to obtain Ij2.6|l . Indeed, by noting that 

X T 2 

Qoidiiv = tAw + -j^oidtw - —9fw 

l<i<3 1<*<3 
„2 



(t,- r —Y 2 t w+ jttoidtw-taw, 



1<!<3 



we see that (12.6(1 follows easily. 
To prove (|2.7|l . we begin with 

(t-r)flS f ft = (t-r)(ifioi-y9 ( ) 



9< 



Since we are in the region r < t/2, which implies that t < ^\t — r\, this yields the bound 

\(t-r)did t w\ < \n 0l d t w\ + \t - r\\d?w\- 
By applying (|2.6|l . we see that (|2.7|l follows. 
Case 2 (\x\ > Here, we introduce the vector field 

(2.9) Q r = ^-Qoi = td r + rd t , 

\x\ 

and begin by noticing that 

(2.10) d t fl r w — d r w — td t d r w + rd 2 w, 

(2.11) d r £l r w — d t w — td 2 w + rd r d t w. 
These immediately yield 

(2.12) r{d t £L r w — d r w) — t(d r Q r w — d t w) = r 2 d 2 w — t 2 d 2 w. 
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In order to prove Q2.5|) . we observe that i|2.12[l is equivalent to 

r(d t Q r w - d r w) - t{d r n r w - d t w) = r 2 U w + t 2 {Aw - d 2 w) - {t 2 - r 2 )Aw. 
Rearranging terms, this is 

(t - r)Aw = — !— r 2 Uw + t 2 (Aw - d 2 w) - r{d t Vl r w - d r w) + t(d r fl r w - d t w) 

t ~\~ T - 

which yields the bound 

\{t - r)Aw\ < rpw\ + t\Aw - d 2 w\ + \d t Cl r w\ + \d r w\ + \d r n r w\ + \d t w\. 

Estimate <|2.5[) follows if we use (|2.8[) and that we are in the region t < 2r. Additionally, 
we require that 

(2.13) \d t n r w\ + \d r n r w\ < \9T a w\, 

\a\<l 

which follows directly from the definition 1)2.9(1 . 

To prove 1)2. we begin similarly by adding and subtracting t 2 Ow to (|2.12|) to get 

r(dtft r w — d r w) — t(d r U r w — dtw) = (r 2 — t 2 )d 2 w + t 2 Uw + t 2 {Aw — d 2 w). 
This, in turn, yields 

(t - r)d?w = — !— \t(d r n r w - d t w) - r(d t n r w - d r w) + t 2 U w + t 2 (Aw - d 2 w) 
t + r l 

Using [)2.8|l and (|2.13|) . we see that H2.6|) follows. 

To prove l|2.7|) . we split the gradient into its radial and angular components, 

V x = -d r t Afl. 

r r A 

Beginning with the radial component, we see that 

(f - r)d t d r = t —^n r d r - t—^dl 
r r 



t — r t — r / 2 I \ 

n r d r - 1 (a — d r — -^n -n). 

r \ r r J 



r 

Using that t < 2r, we see that 

\(t - r)d t d r w\ < \ fl r d r w\ + \(t- r)Aw\ + \d r w\ + |Vfiw| 

< \d r n r w\ + \d t w\ + |(t- r)Aw\ + \Vw\ + \VQ,w\. 

By 1)2.13)1 and (|2.5|l . we see that (|2.7|) holds when W x is replaced by d r . 

To prove (|2.7|) for the angular components of the gradient when t < 2r, the bound 
follows trivially as 

|(t-r)^ Afld t w\ < \d t tlw\, 
which completes the proof. □ 
By arguing as in the proof of Garding's inequality (see 8 , Lemma 3.1), we can obtain 
Corollary 2.5. Let w S C£°(M + xl 3 ). Then, 

(2.14) \\(t-r)dl x T a w(t,.)h< Yl rV(V)l|2+ £ \\{t + r)0^u(t,-)\\ 2 . 

\/3\<\a\+l l/3|<M 
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This estimate will be used to control the last term in (|2.4|) . thus providing our main 
wave equation decay estimate. 

3. Linear decay estimates in R + x M 3 . 

In this section, we gather the linear decay estimates that we shall require for the 
boundaryless cases. There are two main estimates. The first is for the linear wave 
equation and will be used to handle eigenfunctions with zero eigenvalue. The latter is for 
linear Klein-Gordon equations and will be used to estimate the eigenfunctions of nonzero 
modes. 

Our decay estimate for the wave equation is a rather simple consequence of (|2.3|) . 
JZ1, and EFTM . 

Proposition 3.1. Suppose that w S C°°(R x R 3 ) satisfies w{t,x) = 0, t < 2B where B 
is a fixed positive constant. Suppose also that Ow(t, x) — for \x\ > t — B . Then. 

(3.1) (l+i)sup|d t ,*™(t,z)| < J2 Y. sup 2 k \\T a Dw(T, .)||a. 

\a\<2 k r£[2fc-\2*+l]n[2B,i] 

In this section, as in the previous, we are using □ = <9f — Ar3 as we are discussing the 
boundaryless wave equation. 

Proof of Proposition \3.1\ Here, again, we examine the regions |x| < t/2 and |x| > t/2 
separately. In the former case, we may apply (|2 .4|> and Ij2.14|l to sec that 

(1+t) sup \dw(t,x)\< \\r a w'(t,-)\\ 2 + \\(t + r)T a a w (t, .)||a. 

{x:\x\<t/2} |Q , |S1 

By the energy inequality, this is 




If we dyadically decompose [2B, t), it follows immediately that this is 
<EE SU P 2 k \\T a Ow(r, .)||a 

|a|<3 k r 6 [2»-S2*+i]n[2B,t] 

as desired. 

Over | a; | > t/2, the estimate follows similarly. Indeed, we apply l|2.3|l and the energy 
inequality to see that 

(1+t) sup \dw(t,x)\<y2 I \\T a n w (T, -)\\ 2 dT, 

{x:\x\>t/2} \ a \<2^ 2B 

and use a dyadic decomposition in the temporal variable, as above, to see l|3.1|) . □ 

The next estimate will be used to bound the eigenfunctions with non-vanishing eigen- 
value. This non-vanishing eigenvalue will serve as a mass term after we expand, and 
thus, we are left with finding decay estimates for Klein- Gordon equations. As in |ll)|. we 
will rely on the following estimate which is essentially from Hormander [3j (Proposition 
7.3.6). 
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Proposition 3.2. Suppose that w S C°°(R x R 3 ) satisfies w(t,x) — 0, t < 2B where B 
is a fixed positive constant. Suppose also that (□ + fi 2 )w(t, x) — for \x\ > t — B. Then, 

(3.2) (1 +t) 3 / 2 sup \w{t,x)\ < V y sup 2 k \\T a {D + ^ 2 )w{t, -)\\ 2 

M<5 ft rep«.-i,2*+i]n[2B,t] 

where the implicit constant is independent of \i > 1 . 

To prove this, one needs only modify the proof in j2j by using the proper version of 
Lemma 7.3.4 of [3\ This states that if v" + [i 2 v = h in [a, b] C R, then 

1 /' h 

sup < \v(a)\ + \v'(a)\ + - / \h(p)\ dp. 

a<p<b H Ja 

4. Linear decay estimates for waveguides. 

In this section, we show that the estimates of the previous section can be adapted to 
yield estimates in waveguides. In the process, we lose some additional regularity which 
depends on the dimension d of Vl C M. d . 

Before we proceed, we review some basics from spectral theory and elliptic regularity 
theory. We refer the interested reader to, e.g., the texts of Taylor Jl] and Gilbarg and 
Trudinger for more thorough treatments. 

Here, we allow Aq to denote either the Dirichlet Laplacian with boundary conditions 

(4.1) h\ an = 0, 

or the Neumann Laplacian where the boundary conditions are 

(4.2) cUU = 0, 

where d v denotes the normal derivative. Since £1 is compact with smooth boundary, it is 
known that the spectrum of — An is discrete and nonnegative. Letting < X 2 < A3 < . . . 
denote the eigenvalues (counted with multiplicity), we know that 

< \ 2 < Aj < . . . for the Dirichlet Laplacian, 

and 

= A^ < X 2 , < A3 < . . . for the Neumann Laplacian. 

We shall let Ej : L 2 (Vl) — » L 2 (Q) denote the projection onto the jth eigenspace. Thus, 
for h e L 2 (Q), we have that Ejh is smooth and satisfies 

-AnEjhix) = \ 2 E.jh(x). 

Moreover, we have a Plancherel's theorem 

00 

ll^ll!=(a) = Ell^ /l ll2- 

Since by the Weyl formula, Aj « j l / d , j = 2, 3, . . . , we see that for h € C°°(Ci), 

(43) (l+3) 2/d \\E^ (a) <C\\(I-A a )E s h\\L'(a) 

= C\\E ] (I-A n )h\\ L 2 m , j = 1,2,3,... 

provided that either (|4.1|l or 14.2jl holds. 

As for the elliptic regularity theory, we start with the basic estimate 
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Lemma 4.1. Suppose that h £ C°°(fl) and that either l|4.1|l or (|4.2|) holds. Then, for 
N = 2,3,..., 

(4.4) E KHl^u) < E P Q A^|| L2(f2) + INIi^n)- 

H<jV |a|<Af-2 

By writing — An = □ — 9 t 2 + A R 3 and noticing that d t and preserve the boundary 
conditions, it follows from (|4.4I) and an induction argument that 

Corollary 4.2. Suppose that u(t,x,y) S C°°(R+ xK 3 xfi) and that either u(t,x,y)\ y ^dn 
or d v u(t, x, y)\ V £dn = 0. TTien, /or iV = 2, 4, 6, . . . , 

(4.5) E ll 9 X*>*> 011^(0) < E II^X*'^ -)ll^(n) 

|a|<AT |a|<iV 

+ E ll^y^^OIU^O). 
|a|<JV-2 

and for N = 2,3,4,..., 

(4.6) E K«fr*, -)Mn)< E Oilmen) 

q|<JV M + I/3|<AT 

I/3|<1 

+ E H^,y au (*'^ -)IU-(O)- 

|a|<iV-2 

The following lemma shows how our equation behaves under the eigenfunction expan- 
sion. 

Proposition 4.3. Assume that u(t,x,y) € C°°(R+xR 3 xf2) satisfies either u(t,x,y)\ y& gn — 
or d u u(t, x, y)\ y edfi = 0. Then, it follows that 

EjUuit, x, y) = (df - A M 3 + \ 2 3 )E 3 u{t, x, y) 

and 

Y^\\(d 2 t - ^ + \])E 3 u{t,x, 01li*(n) = IXM, Oilmen)- 

The last of our basic elliptic regularity estimates follows from Sobolev's lemma and 
gill- 
Proposition 4.4. Assume that u(t,x,y) € C°°(R+xR 3 xf2) and that either u(t, x,y)\ y£ gn = 
or d u u(t, x, y)\ v edn = 0. TTien, 
(4.7) 

|a?«(*,^l/)l< E ll$>(*>*>-)lli»(n)+ E 11^°^'^ -)ll^(n)- 

H<|/3| + (d+4)/2 \a\<\p\+d/2 

We are now ready to present our main decay estimates. We let 

(4.8) {f } - {r} n {d y } = {d t , d x ,n jk ,d y : o < j < k < 3} 

denote the full set of "admissible" vector fields. 
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Proposition 4.5. Fix B and suppose that u £ C°°(R_|_ x R 3 x Ct) satisfies u(t, x, y) = 
for t < 2B, Du(t, x, y) = for \x\ > t — B, and d v u{t, x, y)\ y <Edtt = 0. Then, 

(4.9) (l + t)\fPd t Mt,x,y)\< E E SU P 2 fe ||f«Du(r, .)||2 

M<|/3| + (5d+16)/2 k re[2"-K2^]n[2B,t] 

+ (i+«) e iif Q n^, oib. 

M<|/3| + (5d+6)/2 

Here, as above, d is the dimension of CI. 

Proof of Proposition \4-5\ The proof of this proposition follows that of Proposition 3.5 
of JOI very closely. Indeed, the only major modification is the use of to estimate 
those terms in our expansion with vanishing eigenvalues. 

By noting that T preserves the Neumann boundary condition and commutes with □, 
it suffices to take T^ 3 = d@. Using (|4.7|) and Plancherel's theorem, we have that 

(4.10) _^ 

\d$d t , x u(t,x,y)\ 2 < E ll$!A*«(*»*. Oilmen) 

H<|/3| + (d+4)/2 

+ E K x ^u(t,x,-)\\ 2 L2{Q) 

H<|/3| + (d+2)/2 

00 

£ E .)||i 3( n) 

|a|<|/3| + (d+4)/2 j=l 

+ E \\®t,x,y au (ti ' )lll 2 (K3 x n)- 

M<|/3| + (d+6)/2 

In the last step, we have also applied Sobolev's lemma. 

In order to gain further control over the first term in the right, we apply (|4.3|l to see 
that 

(l+j) 2/d E \\Ejdl x d t , x u{t,x, .)||z> ( n) 

|a|<|/3|+(d+4)/2 

< E ||^(J-AoK x 9 t , xU (t,x, 011^(0) 

M<|/3| + (d+4)/2 

i$ E \\ E J d ?,xdt,xU(t,X, -)\\ L 2 m 

H<|/3| + (d+8)/2 

+ E ll^ x Du(t,z, -)IU=(n). 

H<|/3| + (d+6)/2 
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In the last step, we have simply written — Aq = □ — &f + A. If we apply this estimate 
recursively, it follows that 

E \\ E i d t,x d t,xu(t,x, -)IU=(n) 

\a\<\P\+(d+4)/2 

~ JlT ,)2 E \\ E 3 d txdt,xU{t,X, -)\\ L 2 (n) 

1 J) |a|<|/3|+(5<£+4)/2 

+ E WEjd^Uu&X, 011^(0). 

H<|,3| + (5d+2)/2 

By Plancherel's theorem and Sobolcv's lemma on R 3 , this implies that 
(4.11) 

OO 

E Ell^AW*^. Oilmen) < E \\d^ x au(t,-^ 



»|<|/3| + (d+4)/2j = l |q|<|/3| + (5£;+6)/2 



E ttx E Pi^xft,*«(*.«. Oil i 2 (o)- 

4=1 1 7J |a|<|/3|+(5d+4)/2 



Combining l|4.10|l and (|4.11l) , it follows that 
(4.12) (l + t)X^(M,2/)| 2 <(l + 2 E KMt, Olli»(R3xn) 

a|<l/3| + (5<i+6)/2 

+ (l + i) 2 E llfc,x£i8JX*,s,0ll 

a|<|/3| + (5d+4)/2 
oo 

+ ( 1 + < ) 2 E7TT7F E ll^«(t,*,OII 

J=2 1 J ' 1 |a|<|/3| + (5d+6)/2 

Here, we apply l|3.1fl to the second term on the right side, which yields 

M<l/3| + (5d+4)/2 

^ E E SU P 2 fc ||r Q (5 t 2 -A)£; lM (r, .)IU»(R-xn) 

M<|/3| + (5d+8)/2 fc re[2*-i,2fc+i]n[2B,t] 
M<|/3| + (5d+8)/2 k re[2«-\2^]n[2B,t] 

The last inequality follows from Proposition ^. 31 
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We can similarly apply i|3.2[l to each of the summands in the last term of (|4.12|) , where 
Xj > c > for each j — 2,3, ... . This yields 

(1 + t) ^ ||^ x tt(*.^-)IU»(n) 

ct l<l/3| + (5d+6)/2 

^ E E SU P 2 k \\T a {d 2 t -A + X))E jU {r, -)IU»(R3 X n) 

|a|<|/3|+(5d+16)/2 ft r6p»- 1 ,2*+i]npB ) t] 

~ E E SU P 2 fc ||r«n W (r, -)IU»(R»xn). 

|a|<|/3| + (5d+16)/2 fc re[2*-i,2»+i]n[2B,t] 

Here, again, the last inequality is a result of Proposition 14.31 

Combining these last two bounds with (|4.12l) completes the proof. □ 

For the case of Dirichlet boundary conditions, a similar bound can be achieved. 

Proposition 4.6. Fix B and suppose that u € C°°(R+ x R 3 x Cl) satisfies u(t, x, y) = 
for t < 2B, Du(t, x, y) = for \x\ > t — B, and u(t, x, y)\ V £dn — 0. Then, 

(4.13) (l + t)\fPu(t,x,y)\< V E SU P 2 k \\f a n u (T, .)||a 

|a|<|/3|+(5d+14)/2 fc re[2*-i,2^+i]n[2B,i] 

+ (l + t) e l|f Q D^, .)||2. 

|a|<|,3| + (5d+5)/2 

This is essentially Proposition 3.5 from ^U|. In fact, as there is no zero mode, an 
estimate with stronger decay is available. In our proof of long time existence, however, 
we would not be able to handle the added power of (1 + t) 1 / 2 that would be necessary 
in the last term of (|4.13|) . The modification required to obtain this weakened version is 
straightforward, and as the proof closely resembles that of the previous proposition, it 
will be omitted. 

Using this, we can obtain the necessary decay for d y u when Q is one-dimensional. 
Indeed, we simply notice that when O is one-dimensional, v = d y u = d v u satisfies the 
Dirichlet wave equation Dv — d y Ou. Thus, we may apply l|4.13|) to obtain the following. 

Corollary 4.7. Fix B and suppose that u E C°°(R+ xl 3 x [a, b]) for two fixed constants 
a < b. Suppose further that u satisfies u(t,x,y) = for t < 2B, Ou(t,x,y) = for 
\x\ > t — B , and d y u{t, x, a) — d y u(t, x, b) = 0. Then, 

(4.14) (l + t)\t^d y u(t,x,y)\< V V sup 2 k \\f a Du(r, 

M<|/3| + 10 k r6[2*-l,2*+l]n[2B,t] 

+ (i + t) l|r Q n«(t, -)ll 2 . 

M<l/3|+6 

5. Energy estimates. 

In this section, we briefly describe the energy estimates for the variable coefficient 
wave equation which we shall require. For 7 jfc S C°° satisfying 

3+d 3+d 

(5-1) £ ^ ^ X / 2 ' E ll^lUji- ( [0,T]xK*xn) < Co, 

j,k— i,j : k=0 
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we study 

< ac, j/) = o, ye on, 

w(t,x,y)=0, t<2B. 
In order to establish energy estimates with Neumann boundary conditions, we are re- 
quired to assume a natural nonlinear compatibility condition on the 7 jfc . This says that 

3+d 

E l jk (t,x,y)^6 k = 0, if y edn, 
(°- 2 > j,k=0 

0=(p,...,0,Vi(y),...Mv)), and £-0 = 0. 
Under this assumption, if u vanishes for large \x\, the standard proof of energy inequalities 
(see, e.g., Proposition 2.1 in Chapter 1 of ^3]) yields 

(5-3) \\dt, x , y u(t, .)|| 2 < / \\F(s, -)\\ 2 ds. 



In the sequel, we shall require a version of l|5.3|l where u is replaced by T^u. Since T 
preserves the boundary condition, it follows from 1)4. 6fl that 

E \\r a d$d tiX , y u(t, -)\\ 2 < E ||r^ )J/ ^ ) -)l|2+ £ l|f Q n^, 

|a| + |/3|<JV |a|<AT |a|<iV-l 

Thus, since 

3+d 3+d 3+d 

ar a u+ E l jk djd k T a u = T a F + E [7 jfe ,r a ]^a fc w+ £ 7 ^[<9 J <9 fc ,r c >, 
j,fc=0 j,fc=o j,fc=o 

it follows easily that for N = 0, 1, 2, . . . 

(5.4) £ ||f«cW^, .)|| 2 < E [ \\T a F(s,-)hds 

\a\<N \a\<N 

3+d „j 3+d „f 

+ E E / l|[7 J VW feU ( s ,.)IM*+ E E / h^W^OlM* 

\a\<N J,k=0 J ° \ a \<N ],k=0 J ° 

+ e nf' 3 ^, -)ii 2 . 

\0\<N-1 

6. Almost global existence. 

In this section, we prove Theorem O and Theorem O We will prove the former 
explicitly and will mention in the course of this proof the modifications that would yield 
the latter. In order to simplify the exposition, we will take d = 1 in Theorem ll.il also. 
The general case only requires additional occurrences of the vector fields (the number 
depending on d) in the definitions l|6.4|l and l|6.10|l below. 

In order to apply our decay estimates, we shift the time variable by 2B and examine 
initial conditions at time t = 2B. Using local existence theory (see, e.g., |S]) if e > is 
sufficiently small and N is sufficiently large in (|1.4f> . it is known that there is a solution 
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to on an interval [2B, 2B + 1]. We will use this to reduce to the case of vanishing 
data (as was similarly done in Fixing r\ S C°°(IR) with rj(t) = 1 for t < 2B + i and 
r)(t) = for t > 2B + 1, we set 

u (t,x,y) = rj(t)u(t,x,y) 

and note that this solves 

□tt = 7]Q(dt iX u,df x u) + [0,rj\u 

with 

d l ,u (t,x 1 y)\ yedn = 0. 
We also note that the local existence results imply that 

(6.1) sup ^ ||f Q u (i, 0lk 2 (K 3 xn) < C Q £. 

1 \<A<N 

We conclude that u solves for 2B < t < T if and only if w — u — uq solves 

(aw = (1 -ri)Q{dt, x (u + w),dl x (u + w)) - P,t]}(uq + w), 

(6.2) ld u w(t,x,v) = 0, yedn, 
[w(t,x,y) = 0, t<2B. 

We solve (|6.2fl on [0, T e ) using an iteration. We set Wq = and define recursively 
to solve 



(6.3) 



Ow k = (1 - rj)Q(dt,x{ua + w k -i),dl x (u + w k )) - [a,r]](u a + Wk), 
d v Wk{t,x, y) = 0, y € dfl, 
w k (t,x,y)=0, t<2B, 



for k= 1,2,3,.... 
We define 

, x w k (t,x,y) 

\a\<30 l/3|<16 



2B<t<T 



Here, in the proof of Theorem 11.21 we need to, of course, allow d y Wk in the last term. 
We will label the two terms in (|6.4(l If, (T) and Ilk (T) respectively. 

Our first goal is to inductively show that for e sufficiently small, we have 
(6.5) M k {T e ) < 4Cie. 

Here, C\ is a uniform constant that is, say, 10 times greater than (1 + Co) times the 
implicit constants appearing in (|4.9|l and l|5.4|l . When k = 1, l|tj.5[l follows from (|1.4|l . the 
well-known local estimates, (|5.4|) . (|4.9|l . and Gronwall's inequality. 
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Assuming that i|6.5[l holds for k replaced by k — 1, we apply l|5.4|) to see that 
(6.6) 



h(T £ )<c f E 53 II^Muo + «*-!)(*, -)IU \\r a d t Ju + w k )(t, -)\\ 2 dt 

|/3|<15 |a|<30 

+ C f ' J2 \\r d t , x (u o + w k ){t, -)||oc \\^ a dtAuo + w k ^)(t, .)\\ 2 dt 



|/8|<15 M<30 



+ c f E l|r"ft,x(«o + «;fc-i)(t, -)IU 51 l|f a <M u o + w fe _i)(t, .)||adt 

l/3|<15 |a|<30 

+ C sup ^ ||f^ ja! (uo + ti>fc-i)(t, -)IU 5Z l|f°fle,*(«o +«>*)(*. - ) || 2 



- - e |j8|<15 M<30 



•C sup 53 \\f d t , x (u o + w k ){t, -)||oo 5Z l|f Q a M K+^-i)(i, 



Ui: |/3|<15 |a|<30 



C sup 53 ||f^, s (uo + rc fc _i)(t, Olloo 51 ||f a 9 M K + ^-i)(i, OII2 



-- s |(9|<15 M<30 



c l£ + 53 / ||r«a^ fc (t, - )|| 2 rfi. 

, ,^ Qn J2B + (l/2) 



|a|<30 

Here, we have set 

3 

1=0 

Notice that since ^ k vanishes if either j = n + lork = n+l (recall that we are looking 
at d = 1), (|5.2[) is trivially satisfied. For the proof of Theorem 11.21 we instead need to 
refer to II. 8|) to see that l|5.2l) holds. 

By using (|6.4|) . |JOJ, the inductive hypothesis, and Gronwall's inequality, we see that 
(6.7) 

h{T e ) < Cye + Ce 2 + CeM k ^(T s ) + CeM k (T £ ) 

+ C[{M k . x {T e )f + M k . x {T e )M k {T e )](l + [ ° J—dt) 

K J2B t + t / 

< Cie + Ce 2 + \QCC 2 e 2 {\ + log(l + T e )) + CeM k {T e ) 

+ 4CC l£ (l + log(l + T s ))M k {T e ). 



In order to bound II k (T e ), we apply (|4.9|) (and notice that the second term in the 
right side of 14. 9|) is controlled by the preceding term). In the proof of Theorem 1 1.21 we 
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would need to also use Q4.14p . This yields 

(6.8) II k (T e ) 

<V sup V V ||f^ t)X ( U o+Wfc-i)(T, -)lloo Y, \\r a dtAuo+Wk)(T, .)||a 

. r 6 [2i-i,iy+i]n|M,r.] |/3| < 13 |Q| < 26 

+V sup 2* V Hf^xCuo+wO^OIloo \\r a d t , x (u +w k ^)(T, -)\\ 2 

j re[2i-\2i+i]n[2B,T e ] |/3| < 13 |q| < 26 

+V sup y V ||f^ l!e («o+«; fc _i)(T, .)IU J2 \\f a d t , x (u +w k ^)(r, -)\\ 2 

j r6pi-i,2i+i]n[2B,r.] |/3| < 13 |Q| < 26 

+ c ie + c7 sup V ||f a a tWfe (t, -)l|2. 

2B+(l/2)<t<2B+l | a |< 26 

Here, we use Q6.4fl. (|6.1|) . and l|6.7|l to see that 

J/ fc (T e ) < 2C l£ + Ce 2 + CeAf fc _i(T e ) + CeM k {T e ) 

+ C(l + log(2 + T e )) [Mfe(T £ )M fc _i(T s ) + (M fc _i(T e )) a ]. 

In addition to that occurring in (|6.7|) . the log(2 + T £ ) term appears due to the simple 
observation that we are summing over 0(log(2 + T e )) choices of j with 2 J_1 < T e . 
Applying the inductive hypothesis, we see that this is 

(6.9) II k (T e ) < 2Cie + Ce 2 + l6CCfe 2 {l + log(2 + T e )) + CeM k (T s ) 

+ 4CCie(l + log(2 + T s ))M k (T £ ). 

If we combine (|6.7|l and l|6.9[) , we see that this yields (|6.5|l provided that e is sufficiently 
small and T e is as in (|1.5|l with k sufficiently small. 

It remains to show that the sequence w k converges for t € [2B,T e ). By setting 



2B<t<T 



(6.10) A k (T) = _sup_( Hf a 9t,xK - wjk_i)(t, Olb 

|a|<30 



+ ^ (1 + i) sup if^dt^wfc - w fc _i)(t, i, 

y)l) 

|/3|<16 ^ 

and arguing as above, it can be shown that 

A k {T e ) < ^A k _ x {T e ), 

This suffices to show convergence and thus finishes the proof. □ 
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